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The switching dynamics of current-biased Bi2Sr2CaCu20g+5 intrinsic Josephson junctions 
(IJJs) was studied to clarify the effect of d-wave superconductivity and the stack structure on 
the switching properties. High quality IJJs were fabricated, and then the temperature depen- 
dence of the switching probability distribution was measured for the first and second switchings. 
Although the standard deviation of the distribution detected for both switchings showed simi- 
lar saturation characteristics with decreasing temperature, the temperature at saturation was 
about 13 times higher for the second switching than for the first switching. The properties 
of the first switching can be explained in terms of a conventional underdamped JJ, that is, 
macroscopic quantum tunneling below the crossover temperature, and thermal activation with 
quality factor of 70±20 above the crossover temperature. In contrast, the relatively higher effec- 
tive temperature for the second switching evaluated from the switching probability distribution 
suggests a dominant thermal activation process under the influence of the self-heating effect 
even at sufficiently low temperature. 
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1. Introduction 

Current-biased Joscphson junctions (JJs) have been 
accepted as an ideal system to explore the dynamics of 
a macroscopic quantum object under controllable en- 
vironments. The dynamics of these junctions can be 
equivalently represented by particle motion in a one- 
dimensional tilted cosine potential. A particle is stati- 
cally confined in the potential when the bias current (J) 
is sufficiently smaller than the critical current (I c ) and 
when the voltage (V) across the junction is zero. When 
the bias current exceeds the switching current (I sw ) , the 
particle escapes from the metastable state either by the 
thermal activation (TA) process or by the macroscopic 
quantum tunneling (MQT) process. Various aspects of 
the switching properties have been examined using metal 
superconductor JJs, and the effect of damping due to dis- 
sipation on the switching dynamics has been clarified in 
numerous studies . 1 ~ 4 Consequently, the quantized en- 
ergy levels of current-biased JJs have recently been con- 
sidered for quantum information device applications . 5 ~ 7 
On the other hand, JJs composed of high-T c supercon- 
ductors (HTSCs), especially Bi 2 Sr2CaCu208+,5 (Bi2212) 
intrinsic Josephson junctions (I JJs), 8 arc promising can- 
didates for various device applications, such as oscilla- 
tors and microwave sensors. The insulating layer of IJJs 
is a naturally-formed atomic layer in the crystal, and 
thus the insulator has high integrity compared with the 
artificial insulators used in metal superconductor junc- 
tions. In addition, the relatively higher plasma frequency 
for IJJs compared with that for metal superconductor 
junctions suggests that the crossover temperature (T*) 
from the classical to the quantum state should be higher 
as well. These features indicate a high expectation for 
applications of IJJs in quantum information processing. 
Actually, the MQT at relatively higher T*s than those 
of metal superconductor junctions was recently reported 
for HTSC JJs. 9-11 However, IJJs have several inherent 
unique characteristics such as the <i-wave superconduc- 
tivity ; 12 < 13 the self-heating effect } 14 ~ 18 the multiple 
stack structure, and the strong electron correlation effect. 
How these characteristics affect the switching dynam- 
ics both in the quantum and classical regimes remains 
unclear. To clarify this, detailed studies on the switch- 
ing properties of IJJs, especially the damping effect, are 
needed. 

Here, the switching dynamics of IJJs was clarified by 
measuring the escape rate using the ramped bias current 
method. One key issue in quantum device applications of 
HTSCs is how the peculiarities of IJJs contribute to the 
damping in the switching dynamics. Here, the switching 
probability distribution P(I) was measured for both the 
first and second switchings in a multiple-branched struc- 
ture in the current- voltage characteristics. The switching 
probability distribution measured for the first switching 
showed the conventional crossover from the TA to the 
MQT. The quality (Q) factor of the first switching was 
then evaluated by applying statistical analysis to the TA 
process at a temperature (T) above the crossover temper- 
ature T\*. The results clearly indicate undcrdamped fea- 
tures of IJJs. The temperature dependence of the switch- 
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ing properties in the MQT regime observed here suggests 
that nodal quasiparticlcs (QPs) have no significant effect 
on the damping in the switching dynamics. In contrast 
to the switching properties of the first switching, those 
of the second switching 19 significantly differed from con- 
ventional features of an underdamped JJs. Although 07 
defined by the standard deviation of P(I) became sat- 
urated for the second switching with decreasing tem- 
perature, the detected crossover temperature T 2 * °f 6.5K 
was significantly higher than conventional metal super- 
conductor cases. We also found the reversal of the rel- 
ative magnitudes of the currents corresponding to the 
first and second switchings, that is, the mean switching 
current for the second switching {I SW 2) was higher than 
that for the first switching (I sw i) when the temperature 
was higher than the crossover temperature (Tcr—1 -IK), 
while I SW 2 was lower than I sw \ when the temperature 
was lower than Tcr- The self- heating effect due to the 
dissipative current inside the junction seriously affected 
the switching properties, and the self-heating assisted TA 
dominated the switching process in the second switching 
below T*. 

2. Model 

In this section, we briefly introduce a model of switch- 
ing dynamics of IJJs by focusing on the difference be- 
tween the IJJs and metal superconductor junctions. The 
JJ is described by the resistively and capacitively shunted 
junction (RCSJ) in which R is the shunt resistance and 
C is the shunt capacitance in the Stewart-Mccumber 
model . 20 Because the IJJs fabricated in this study had 
seven JJs as described in the next section, we assumed 
a serially-connected RCSJ as shown in Fig. 1(a). In the 
model by Koyama et al. 21 each junction is coupled with 
neighboring junctions due to the charging effect induced 
by insufficient shielding of the electric field at the atomic 
size electrodes. If we neglect the effect of coupling, the dy- 
namics of a single junction can be represented by damped 
particle motion in a one-dimensional tilted cosine po- 
tential where the effective mass is given by C($o/2tt) 2 , 
the damping coefficient by 1/ RC , and the magnetic flux 
quantum by <&o - 22 For a constant I smaller than I c , the 
barrier height Uq is given by 

Uq = 2Ej [ \/l — x 2 — x arccos x] , (1) 

where Ej = ^qI c /2tt and x = I/I c . At T higher than 
T*, the escape rate of the particle is dominated by the 
TA process. In this case the escape rate is given by 

r TA = a t ^exp(-^). (2) 

Here, at is the temperature- and damping-dependent 
thermal prefactor and u>o(= ^/27r/ c /$oC(l — x 2 ) 1 ^) is 
the oscillation frequency at the local minima of the po- 
tential , 23 ' 24 The escape rate below T* (= hjj /2-Kk B ) 
is dominated by the MQT process, whose escape rate is 
given by 

/ B \ 1/2 

Vmqt = VeO^o — cxp(-B), (3) 
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where 



= 5^(1 + 0.87/3), 



(4) 



and Q = uj^RC . 25 The experimental escape rate F(I) 
can be written as a function of the experimental P(I) as 
follows, 



r ( i)^iin.^ W 



(5) 



dtAI ir+AiPindi' 1 

where dl/dt is the sweep rate of the bias current. The 
escape rate in the MQT regime is insensitive to the bath 
temperature when the dissipation during the tunneling 
process is negligible, whereas when ohmic (super-ohmic) 
dissipation is non-negligible, the correction to the escape 
rate has T 2 (T 3 ) dependence. 26 Several theories have 
been presented to estimate the effect of nodal QPs of 
d-wave superconductivity on the damping. 27-30 Accord- 
ing to Kawabata et al. , 30 nodal QPs induce superohmic 
dissipation on the MQT process reflecting the QP spec- 
trum of the d-wave superconductors, although the cor- 
rection is quite small. This theory by Kawabata et al., 
however, needs to be confirmed experimentally. Moreover 
it is still unclear how the switching dynamics is modi- 
fied by the coupling of neighboring junctions. Therefore, 
detailed study on the swhitching dynamics of IJJ is re- 
quired. 

3. Experimental 

The IJJs in this study were fabricated as follows. First, 
a single crystal of Bi2212 with a critical temperature of 
90K was grown by using the traveling solvent floating 
zone method. Next, the crystal was cut into 1000x20x5 
/xm size samples by mechanical polishing and cleaving. 
The samples were then glued onto SiC>2 substrates. Four 
Au contact electrodes were deposited onto both ends of 
each sample to achieve an electrical contact to Bi2212. A 
bridge structure was formed on each sample by "necking" 
the center region, and then an IJJ was formed by cut- 
ting the side of the sample by using a focused ion beam 
(FIB). 31 The fabricated IJJ size was 1 x 2 x O.Ol^m 3 , 
and seven JJs were stacked together to form an IJJ. To 
minimize degradation, the IJJ was immediately fixed on 
a sample holder installed in a dilution refrigerator. The 
critical temperature (T c ) of the present IJJ was about 80 
K. Figure 1(b) shows a schematic of the junction config- 
uration and a scanning ion microscope image of the IJJ. 
Based on the size of the IJJ, the estimated capacitance 
was between 74 (relative dielectrical constant e r =5) and 
147 fF (e r =10). All the electrodes were apart more than 
250/im from the junction, and thus the contact resistance 
was less than 10S7 32 . The advantage of this stacked struc- 
ture over the mesa structure 33 is that the effect of the 
heat generated at the contacts is negligible. Actually, no 
dependence of the switching probability distribution on 
such self-heating was detected when the duty ratio of 
the bias signal was decreased from 1/10 (typical value) 
to 1/80 at 52mK. On the other hand, the stacked struc- 
ture is rather sensitive to the heat generation because 
the heat-sinking channel is limited to narrow and poorly 



conducting arms holding the stack. This device struc- 
ture results in relatively high effective temperature at 
the second switching as described below. The quality of 
the junction was confirmed by large hysteresis in the I-V 
characteristics as shown in Fig. 2(a). 

The switching probability distribution was measured 
by using the time domain measurement with the ramped 
bias technique . 1_4 Although the IJJ had four terminal 
configurations, negative signal lines (V- and /_) were di- 
rectly connected to the sample holder as shown in Fig. 
2(b). To reduce the effect of external electric noise, signal 
lines were carefully filtered using low-pass filters (LPF) 
located at room temperature, high attenuation coaxial 
cables, and powder microwave filters (PF) installed in 
the sample holder. The cutoff frequency of the wiring 
was 1 MHz (or 10MHz), and an attenuation greater than 
120dB was achieved at 1GHz. A ramped bias signal was 
generated by using an Agilent 33220A waveform gener- 
ator, and fed to the IJJ through an isolation amplifier 
and a bias resistor Rj. The switching probability distri- 
bution was detected by using a Stanford Research time 
interval counter SR620 in which the time interval be- 
tween the start and stop pulses from the comparators 
was measured. The time interval directly corresponds 
to the switching voltage of the bias signal because the 
bias voltage was swept at a constant rate. The distri- 
bution of the voltage was converted to that of current 
through the I-V relation that was precisely measured 
just before the switching probability distribution mea- 
surements. Primary result for the first switching was 
shown in Ref. 32. For the measurements of second switch- 
ing, we adjusted the sweep rate, the comparator thresh- 
old voltage and the preamplifier gain. The typical sweep 
rate of the bias current was ~ 0.01 lA/s and ~ 0.013A/s 
for the first and second, respectively, and the typical sam- 
pling number to obtain one probability distribution was 
between 5000 and 20000. The threshold voltage of the 
comparator was 2mV(= Vthi) an d 20mV(= VtKi) for 
the first and the second switching, respectively (see Fig. 
2(c)). In determining the sweep rate of the current for the 
second switching, we approximated the I-V characteris- 
tics of the first branch as a linear curve near the switching 
region (typically 18-21/iA). Estimated error due to this 
approximation was less than 1%. The time interval of 
about 0.2~lx 10 _3 [sec] between the first switching and 
the second switching is enough to avoid the dynamical 
effects. We have checked it by changing the sweep rate 
of the bias voltage. All analog circuits were isolated from 
the digital apparatus by optical fiber links to isolate the 
circuits from the digital noise sources and to avoid the 
ground loop. The temperature of the sample Tb at h was 
monitored by using a RuC>2 thermometer tightly fixed to 
a sample holder made of oxygen-free copper. Deviation 
in Tb a th from the actual sample temperature was min- 
imized by fixing the thermometer to the sample holder 
just behind the sample. Further details of the measure- 
ment setup were reported elsewhere. 3 

4. Results 

The I-V characteristics of the IJJs exhibited the 
multiple-branched structure due to successive switching 
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of the stacked JJs in the IJJ from zero to finite voltage 
states as shown in Fig. 2(a). In this study, we focused 
on the first and second switchings with mean switch- 
ing currents of I sw i and I sw %, respectively. The shape of 
P{I) clearly tended to be sharper with decreasing tem- 
perature. 32 This trend saturated around 0.5K, as clearly 
observed in the temperature dependence of 07 shown in 
Fig. 3. Figure 3 also shows a plot of 07 for the second 
switching. When the temperature was higher than the 
crossover temperature (7\* ~ 0.5 A" for the first switching 
and T 2 * ~ 6. 5 A' for the second), both switchings followed 
a T 2 / 3 dependence, indicating a dominant TA process in 
the switching. The saturation at T* in the first switching 
indicates the crossover of switching mechanism from the 
TA to the MQT , 35 Actually, the measured T* is simi- 
lar to previously reported values. 10, 11 On the other hand, 
the significantly higher T 2 than T£ suggests the presence 
of different origins for the second switching. Since eval- 
uated temperature dependence of the critical current for 
first switching 7 cl and for second switching J c2 suggests 
the presence of interaction between the two junctions, we 
will discuss this topic in section 4.2. 

4-1 First switching 

The dynamics of the first switching were analyzed here 
to determine the magnitude of the damping from exper- 
imental results shown in Fig. 3. Similar saturation be- 
havior of 07 due to the MQT and relatively weak dissi- 
pation behavior in Bi2212 IJJs have already reported in 
rcf. 10. Here we mainly focus on two points which have 
not been clarified in the previous work; i) the damping 
effect evaluated from the switching properties above the 
crossover temperature, ii) the influence of nodal QPs of 
the d-wave superconductivity on the switching properties 
from the switching properties below the crossover tem- 
perature with high precision measurements. As a whole, 
we will present a consistent analysis on the switching 
dynamics of the IJJs in terms of the underdamped JJ 
behaviors for the whole temperature range. 

Clarifying the damping effect is crucial for understand- 
ing the dynamics of a J J. According to the Caldeira- 
Leggett formula , 25 the MQT rate is suppressed by the 
damping effect. Thus the magnitude of the damping can 
be estimated from the crossover temperature and the 
saturation value of 07 - 1 However, such analysis is not 
adequate for the present case because the estimation of 
Q factor based on this method is sensitive to the capaci- 
tance, and the capacitance is roughly estimated from the 
size of the junction. Another possible way to estimate the 
Q factor is to evaluate the gradient of I-V characteris- 
tic inside the gap. However this evaluation cannot be 
validated, because the experimental I-V characteristics 
is influenced by the AC Josephson effect in the absence 
of applied field. 35,36 Also the gradient of I-V near zero- 
bias voltage is hard to be evaluated due to the retrapping 
effect. Moreover, the relation between the I-V character- 
istics and the Q factor is unclear when the I-V charac- 
teristics inside the gap show strong non-linearity due to 
the ci-wave superconductivity, which is the condition in 
our study. Therefore, our approach was to extract the Q 
factor from the TA process in which the escape rate de- 



pends on the damping, because the escape rate reflects 
the distribution of particles confined in the cosine po- 
tential. Based on the analyses described in refs. 2,4, the 
escape temperature T esc of the TA for the first switch- 
ing can be calculated by setting a t = 1 in Eq. 2 and 
assuming the cosine potential of Eq. I. 37 Figure 4 shows 
the Tb a th dependence of T esc obtained using the least 
squares method. Below Tj", T esc was constant, reflect- 
ing the temperature-independent switching of the MQT. 
Above 2]*, although T esc was approximately Tb at h, there 
was a non-negligible discrepancy between T^ath and T esc , 
and this discrepancy tended to increase with increasing 
Tbath- This discrepancy was not due to measurement er- 
ror or to heat generation at the contact electrodes. Such 
discrepancy is similar to that observed in high-Q junc- 
tions of metal superconductors reported in ref. 4. The 
most plausible origin is the correction in the escape rate 
due to the deviation in a t from unity when the IJJ is 
in the underdamped regime. The Q factor of the junc- 
tion can be evaluated by assuming the following rela- 
tion, 4, 23,24 



(Vl + (QW1.8(7 ) + 1) 2 

Figure 4(b) shows the obtained Q factor in the tem- 
perature range from 0.6K to 5K. Although the results 
were somewhat scattered due to statistical error, the 
Q factor was approximately 70 ± 20, indicating that 
the IJJ was clearly in the underdamped regime. This 
Q factor is somewhat higher than that obtained by 
an YBa 2 Cu3 07_5 junction evaluated using microwave- 
assisted tunneling, but comparable to that obtained by 
high quality metal superconductor JJs. 3 ' 4,38 

In the following, the dynamics of the first switching is 
shown to be consistent with that of a conventional un- 
derdamped JJ without the influence of the nodal QPs 
of the ci-wave for all temperature range below 5K. First 
of all, the Q factor above the crossover temperature Tj* 
should be an exponentially decaying function of temper- 
ature if the Q factor is originated by the nodal QPs of 
the ci-wave. 39 This is because the amount of QPs con- 
tributing to the damping increases with increasing tem- 
perature. The expected Q factor in such situation is the 
order of 10 5 at 0.5K according to rcf. 39. On the other 
hand, the experimentally detected Q factor shown in Fig. 
4(b) is approximated by by Q=69.3+0.8T. This value is 
high enough to be regarded as underdamped, but both 
the magnitude and the temperature dependence are quite 
different from those expected for the nodal QPs damping. 
Secondly, it is accepted that the crossover temperature 
is the function of damping as shown in ref. 25. The ex- 
perimental result on Tj" of about 0.5K is consistent with 
the underdamped theoretical value of 0.445K. Thirdly, 
although the Q factor below T-j* is not directly evaluated 
in the present experiments, o7=50-52nA are consistent 
with estimated value of 51.1nA (7 c =21.6/zA, C=105fF, 
e r =7.16) obtained from the underdamped MQT the- 
ory. 25 Also the lack of the temperature dependence of 
the a 1 is inconsistent with 07 (T) — 07(0) oc T 3 which is 
the theoretical temperature dependence if the nodal QPs 
of the d-wave is the dominant origin of damping. 30 This 
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fact indicates that the nodal QPs are not the dominant 
origin of the damping in the MQT region. The lack of 
the dominant influence of the nodal QPs of the d-wave 
on the MQT process is consistent with the theoretical 
estimation of ref. 30. As the whole, the dynamics of IJJ 
is consistent with that of an underdamped IJJ, and the 
dominant origin of the residual damping is not due to 
nodal QPs but to extrinsic factors such as microshorts 
in the IJJ or coupling to the environment via the signal 
lines. Therefore, the Q factor can be improved by refin- 
ing the junction fabrication process and by isolating the 
signal lines with proper microwave engineering similar 
to the technique used to improve the Q factor in metal 
superconductor junctions. 6,40 

4-2 Second switching 

In this section the properties of the second switching 
were studied by taking account of the self-heating ef- 
fect and of the stack structure. Although the junctions 
in the IJJs are capacitively coupled with each other, 21 
it is still unclear how this coupling affects the switch- 
ing dynamics of each junction. 22 This is an important is- 
sue because weakly coupled J Js arc promising candidates 
for multi-qubit applications. 40,41 We analyze the switch- 
ing current distribution under the influence of dissipa- 
tive current based on the conventional TA model . 23,24 
The advantages of the present method are simple device 
structure and high sensitivity of the switching current 
distribution on the temperature. 

Three notable features were observed in the result of 
second switching and details are discussed below: 

1. The detected crossover temperature (T 2 *) was 6.5K 
(Fig. 5) which is far higher than theoretically expected 
value of 0.74K. 

2. The reversal of the relative magnitude of the mean 
switching current was detected, that is, I sw \ > I SW 2 
below 1.1K and I sw \ < I SW 2 above 1.1K (Fig. 6). The 
relation of I sw i > Isw2 is inconsistent with the naive 
expectation that each junction in the stack switches to 
the finite voltage state on the order of smaller I c when 
the bias current is gradually increased. 

3. Although it is accepted that the temperature de- 
pendence of I c follow AB formula in the cases of IJJs, 
both of the temperature dependences of I c \ and J c2 in 
Fig. 6 are inconsistent with Ambcgaokar-Baratoff (AB) 
formula. 

For the barrier potential of the second switching, we 
assumed a simple tilted cosine potential as described in 
Eq. 1. Figure 5 shows the temperature dependence of 
T esC 2 for the second switching obtained using the least 
mean squares method. The T esC 2 shows saturation be- 
havior at about T 2 * ~ 6.5A'. Because the calculated T 2 
from the TA to the MQT for the second switching is 

0.74K (=5hJ^k/36k B ), the observed T* is too high 
to expect the transition to the quantum regime. Further- 
more, even above T 2 *, significant discrepancy is evident 
between Tbath and T esC 2- Such large discrepancy is incon- 
sistent with conventional JJ behavior, even assuming un- 
derdamped conventional JJs. However, these anomalies 



are understood in terms of the self-heating effect, which 
is due to the poor thermal conductivity of Bi2212. 14-18 
When all the junctions are in the zero-voltage state, no 
dissipativc current exists in the IJJ, and thus the effec- 
tive electron temperature in the stack should be equal 
to the Tbath- Whereas the influence of the dissipativc 
current cannot be neglected when the first junction has 
been transited to the finite voltage states. The most 
plausible explanation is that the measured T esc2 directly 
corresponds to the effective electron temperature in the 
stack and the switching events are induced by the en- 
hanced thermal fluctuation (referred to as self-heating 
explanation). 42 Along this explanation, the evaluation 
of switching distribution can be regarded as an in-situ 
measurement of the effective electron temperature in the 
stack. 17, 18 Actually, similar discrepancy between T esC 2 
and Tbath has been reported for the stacked structure 
as well as mesa structures in which QPs were injected 
into the JJ from a metal electrode placed near the junc- 
tion. 37 The uniqueness of the present work is that the 
sensor used to detect the effective temperature (secondly 
switched junction) locates quite close to the heating spot 
(firstly switched junction), and then the influences of 
noncquilibrium effect cannot be neglected in addition to 
the diffusive joule heating. 

In order to compare the present results with those 
in previous works, we evaluate the thermal resistance 
R t h = AT/P where AT and P are the change of tem- 
perature and the injected power, respectively. The value 
of Rth in our measurement is given by 2.3 xl0 7 (K/W) 
(AT -8.5K, and P ~ 17 mVx 21.6/iA), which is far 
higher than those detected in mesa . 18 Also the AT es- 
timated based on the equation presented in ref. 16 tak- 
ing account of mesa structure gives 0.36K by assuming 
K a b at <1K is lW/Km and K a b/n c ~10 , 43 which is far 
smaller than the present results. These serious discrep- 
ancies reflect the high sensitivity of the stacked junc- 
tion structure for the self-heating as compared to that of 
mesa. Whereas in the case of the stacked structure simi- 
lar to that used in the present work, 17 the relatively high 
thermal resistance of ~2.8x 10 5 (K/W) with AT ~31K 
(110/iW) was obtained. Even taking account of the facts 
that k c at 35K is about 3 times larger than k c at 8.5K, 
and that the distance between the sensor and the heat 
spot in ref. 17 is estimated more than ten times as far as 
that in our case, the evaluated thermal resistance in our 
experiments is effectively several times larger than that 
obtained in ref. 17. However, such difference is not unrea- 
sonable because the device used in ref. 17 implemented 
two side bars and they worked as additional heat sinking 
paths. Another important factor for the relatively en- 
hanced self-heating effect in our measurement may stem 
from the position of the " temperature sensor" . In the 
present measurement, the temperature was monitored by 
the secondly switched junction that locates quite close 
(1.5-9nm) to the firstly switched junction. Therefore, in 
addition to the influence of the simple joule heating, 
serious influences by the ballistic phonon and by the 
noncquilibrium quasiparticle distribution (NQD) due to 
quasiparticle injection at firstly switched junction cannot 
be neglected. This is because the distance of those two 
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junctions was less than the c-axis quasiparticle relaxation 
length (about 16nm estimated in ref. 44) as well as the 
mean free path of phonons. Although the influence of the 
NQD on the switching current distribution is not clear 
at present, it can be approximated as elevated electron 
temperature similarly to the usual treatment of noncqui- 
librium superconductivity . 45 Therefore, we ascribe that 
the origin of the enhanced AT in the second switching 
to the self-heating effects, that is, the sum of the simple 
Joule heating, the ballistic phonon and NQD. 46 

The TA process analysis treats the mean switching 
current (I S wii Isw2) and the critical current (I c i, I C 2) as 
independent parameters. Figure 6 summarizes the analy- 
sis results for the IJJs studied here. At the first switching, 
the difference between I c \ and I sw \ was less than 1/iA be- 
low Tj", and the difference rapidly increased when Tf, Qt ^ 
was increased from Tj*. In contrast, at the second switch- 
ing, the difference between I C 2 and I SW 2 was almost at 
6/iA for all Tt, at h below 5K because of the increased effec- 
tive electron temperature due to the heating effect. The 
difference in the temperature dependence between I sw \ 
and I SW 2 causes the reversal of the relative magnitudes 
of currents between the first and second switchings at 
about I .IK (=Tcr)- The fact that I sw i > I SW 2 has also 
been reported in ref. 10. 

Here we discuss the anomaly of the temperature de- 
pendences of I c i and I C 2 in Fig. 6. It is accepted 
that the temperature dependence of I c follows AB for- 
mula as suggested by Krasnov et al. 37 Based on the 
AB formula, the critical current keeps constant value 
(0.9999 < I C (T)/I C (T = 0) < 1.0) in the tempera- 
ture range < T < T c /5, therefore I c should be con- 
stant in the present experimental temperature range 
shown in Fig. 6. However, I C 2 in Fig. 6 is gradu- 
ally decreasing as the temperature is increased. The 
amount of change is / c2 (T=5K)// c2 (T=60mK)=0.969, 
which is far larger than that expected from AB for- 
mula. Furthermore, I c \ in Fig. 6 is increasing as the 
temperature is increased with the amount of change of 
7 c i(T=5K)// c i(T=60mK) = 1.016. The fact that the I c 
is an increasing function of temperature is anomalous 
because it is inconsistent with conventional theories for 
weak-link Josephson junction 47 as well as AB formula. 
One possible explanation for these anomalies is that I c 
of the IJJ in the stack is modified from that of isolated 
IJJ due to the coupling effect, 21 thus it could have differ- 
ent temperature dependence from that of isolated IJJ. 22 
However we don't have detailed model to explain the 
mechanism of the anomalous temperature dependence of 
I c in the stacked junction at present. 

Apart from above-mentioned anomalies, the dynam- 
ics of the second switching is rather consistent with the 
conventional TA model of an isolated JJ. As shown in 
Fig. 7, P(I) in the second switching consistently follows 
the theoretical curves of the TA model, and no additional 
structures, such as peaks and dips, were found within the 
present resolution. This consistency suggests that the os- 
cillation of the macroscopic phase (AC Josephson effect) 
at the firstly switched junction does not seriously modi- 
fies the second switching properties. On the other hand, 
the influence of the self-heating effect is serious for mul- 



tiqubit application of the junction stack because the self- 
heating might destroy quantum coherence of the others. 
It is important to check whether this problem is peculiar 
to the present stacked structure or not. The Joule heat- 
ing may be reduced in the case of mesa structure because 
of the presence of the heat sinking path at the topmost 
Au layer. Also the miniaturization of the junction size 
would be effective in order to reduce the influence of the 
Joule heating. If we adopt the similar size dependence for 
the present junction to that of mesa, 16 the junction size 
of 47nm gives the heat generation of 0.5K by assuming 
that the current density is size independent. However, it 
is not clear at present whether the self-heating can be 
reduced by the miniaturization of the junction size or 
the implementation of the additional heat-sinking path, 
because the non-equilibrium effect cannot be avoided by 
these modifications. Further studies are needed in order 
to realize the multiqubit in a single IJJs. 

5. Summary 

In this study, high quality Bi2Sr2CaCu20s+5 intrinsic 
Josephson junctions (IJJs) with large hysteresis in the I- 
V characteristics were fabricated. The escape rate mea- 
surements for the IJJs clarified the switching dynamics of 
the IJJs as follows. At the first switching, both the escape 
temperature and the width of the switching probability 
distribution saturated below 0.5K, indicating the transi- 
tion of the switching mechanism from the thermal activa- 
tion (TA) to the macroscopic quantum tunneling (MQT). 
The high quality (Q) factor (70±20) obtained for the first 
switching reveals underdamped features, despite possible 
damping due to the nodal quasiparticles of the d-wavc 
and the stack structure. Because the dominant origin of 
the damping is probably not due to the inherent proper- 
ties of the IJJs, higher Q factor might be achieved by 
improving the connection of the signal lines by using 
proper microwave engineering. At the second switching, 
however, the escape temperature is almost temperature 
independent below 6.5K, indicating significant influence 
of the self-heating effect. The switching of the first junc- 
tion into a finite voltage state increased the effective elec- 
tron temperature, and the escape temperature directly 
corresponds to the effective electron temperature in the 
stack. Although the self-heating explanation reasonably 
describes the enhanced fluctuation in second switching, 
the anomalous behaviors in the temperature dependence 
of critical current remain unsolved. The heating effect 
might destroy the quantum coherence in multiqubit of 
IJJs, however, this problem will be overcome if the IJJ 
with junction size of several tens of nanometers would 
be fabricated. In conclusion, the results from this study 
strongly indicate that high-T c superconductor junctions 
can be used in quantum information device applications. 
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Fig. 2. (color online), (a) Schematic and photograph of the I-V 
characteristics of the stacked IJJ at 55mK. Large hysteresis con- 
firms the high quality of the IJJ. The first and second switchings 
shown here respectively correspond to the transition from the 
zero to the first branch and from the first to second branch. I sw ± 
and I sw 2 indicate the switching current for the first and sec- 
ond switchings, respectively. The relation I sw \ > I S w2 is clearly 
evident in the I-V characteristics (see text), (b) Schematic of 
the electronics used to measure the first and second switchings. 
Bias current was introduced from the bias circuit to the junction 
through a bias resister (Ri) and filters. Both the current bias line 
(1+) and voltage monitor line (V+) were carefully filtered by a 
low-pass filter (LPF) and powder filter (PF). Other connections 
(I- , V— ) were connected to the sample holder, (c) Schematic of 
time-oscillograms used in the present measurements. 



Fig. 1. (a) Equivalent circuit model of the stacked IJJ structure. 
The stacked IJJ includes seven JJs in a stack, and the RCSJ 
model represents each junction, (b) Schematic and scanning ion 
micrograph image of the stacked IJJ structure. 
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Fig. 3. (color online). Temperature dependence of 07 defined by 
the standard deviation of the switching probability distribution 
P(I) for the first switching (red symbols) and second switching 
(blue symbols). Solid line represents a guideline for T 2 / 3 depen- 
dence, and dotted line represents the theoretical value (51.1nA) 
for the MQT without damping. 




Fig. 4. (color online). Bath temperature dependence of the es- 
cape temperature for the first switching (T e3C i) based on the 
TA model with at=l. Solid line corresponds to T esc = Tj, ai ;,. 
Inset (b) shows bath temperature dependence of the Q factor 
for the first switching based on the TA model with correction 
of the thermal prefactor taken into account (see text). Dotted 
line represents fitting to the distribution (Q=69.3+0.STi >at f l , in- 
dicating that the Q factor was relatively independent of bath 
temperature. Inset (c) shows switching probability distribution 
P(I) at T bath = lOmK and the theoretical fitting (J c =21.637/jA, 
C=105fF, Q=70) based on Eq. (3). 
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Fig. 5. (color online). Bath temperature dependence of the es- 
cape temperature for the second switching (T esc 2). Solid line 
corresponds to T esc 2 = T^ at ] x , and the dashed line corresponds 
to the saturation temperature (T esc 2 — 8.5K) of the escape tem- 



Fig. 6. (color online). Bath temperature dependence of the mean 
switching current I sw \ and I sw 2 extracted from the experimental 
data and the critical current I c i and I c 2 from the fit to the TA 
model for the first switching (red symbols) and second switching 
(blue symbols). The crossover of the magnitude of I c \ and I c 2 
occurs at Tcr (— 1.1K ). 
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Fig. 7. (color online). Switching probability distribution P(I) for 
the second switching in the temperature range between 0.55K 
and 9.57K. Symbols represent experimental data, and solid lines 
represent the theoretical fitting based on Eq. (2). 



